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Abstract. We perform a detailed study of inverse Compton (IC) emission for a fireball undergoing external shock (ES) in 
either a uniform or a wind-like interstellar medium, and assess the relative importance of IC and synchrotron emissions. We 
determine the primary model parameters driving the IC to synchrotron emission ratio in the case of a short duration central 
engine. We then investigate the case of ES by a long duration central engine, or delayed external shock (DES), a model that 
can account for some of the flares observed in gamma-ray bursts (GRB) X-ray light curves at late times. We present model 
predictions, in particular in terms of GeV vs X-ray behavior, and compare them with other models proposed to explain the 
origin of flares. We find that if most of the emission occurs when the fireball is in the fast cooling regime, then a substantial 
GeV emission is expected both for a short (standard ES) and a long (DES) duration central engine activity. In particular, in the 
context of standard ES we are able to account for the delayed emission observed in GRB940217. In the case of DES, we find 
that IC scattering of X-ray flare photons can produce high energy flares in the GeV band, which can be detected by the Large 
Area Telescope aboard GLAST. The detectability of high energy flares improves with the burst kinetic energy E: about 30 % of 
Swift GRBs showing flares in their X-ray light curve have sufficiently large kinetic energy so that the expected high flares can 
be detected by GLAST. One important prediction of the DES model is the simultaneity between low and high energy flares. 
To test this simultaneity, the peak energies of both flares need to fall below or within the observational bands. We predict that 
X-ray flares with peak energy of ~ 10 eV produce high energy flares with peak energy of around 100 MeV-GeV. Observations 
by Swift and GLAST then, can test the predicted simultaneity, thus helping to discriminate between different models. 
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1. Introduction 

One of the key areas that would clarify the nature of gamma- 
ray burst (GRB) progenitors is to understand the mechanism(s) 
acting during the GRB prompt-to-afterglow transition phase. 
This phase takes place from hundreds to thousands of sec- 
onds after the burst, and is characterized by a variety of tem- 
poral and spectral behavior due to the contribution of both 
prompt and afterglow emissions. Particularly interesting is the 
appearance of X-ray flares. The comprehension of this phe- 
nomenon can give important information on the physics of 
the central engine. Swift observations showed that X-ray flares 
are present in a large fraction of GRBs (about one half of the 
Swift GRB sample) and occur both in long and short bursts 
dO'Brien et al. 1 2006). X-ray flares exhibit a different behavior 
in terms of their spectral properties: most of them present hard- 
to-so ft spectral evolution resembling that of the prom pt emis- 



Piro et al. I d2005l) : iGalli & Piro I (I2006I) ). Several models have 
been proposed in the literature to explain the origin of GRB 
flaring activity. The most important difference between these 
models is that some of them do not require a long dura- 
tion central engine activity (group 1), while other require a 
long duration and/or re-activation of the central engine (group 
2). Models which belong to group 1 are l ate internal shocks 
(LIS) from a sho rt duration central en gine (|Zhang et a~ 2006; 



1998 



mpt 

sion dBurrows et al. 1 120051 : iButler & Kocevski 1 120071) . while 



others seem to exhibit a softer and not evolving spectrum con- 
sistent with that of the subsequent afterglow emission (e.g. 
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Wu et al. 200a) . refreshed shocks dRees & Meszaros 
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Kumar & Piran 2000b; Guetta et al. 2007), two components 

jet (Meszaros & Rees 2001; ILipunov et al. 11200 ll) . patchy jet 
dKumar & Piran 2000al) . forward shock (FS)-reverse shock 
(RS) dFan& Weill2005h. external shock (E S) with a clumpy 
medium ( Dermer et al. 1 2000: De rmer 120071) and delayed mag- 
netic dissi pation in strongly magnetized ejecta caused by exter- 
nal shock (Gia nnios 1 l2006h . 

In the scenario of LIS from a short duration central en- 
gine, in addition to the shells producing the prompt-gamma 
ray emission through internal shocks, a tail of slow (small 
Lorentz factor) shells is emitted by the central engine. These 
slower shells can collide (internal dissipation) at a later time 
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thus producing (X-ray) flares. In this case a late prompt emis- 
sion is produced, while there is no external shock emission 
because the medium has been previ ously cleaned up by the 
faster shells (ILazzati & Perna 1 12006). In this context, in or- 
der to explain a delayed X-ray emission appearing hundreds 
to thousands of seconds after the burst, a small difference in 
the Lorentz factor F of the colliding shells would be necessary, 
i.e. Ar << T. However a small contrast between shells can- 
not explain the f ast temporal v a riabili ty typically observed in 



X-ray flares (see iKrimm et al. I (|2007) for details). Refreshed 



shocks occurs when the central engine releases its energy with 
a variety of Lorentz factors, thus the faster part of the outflow 
decelerates earlier and the slower part can catch up with it at 
later times (when it has decelerated due to the interaction with 
the external medium), injecting energy into the blast wave and 
producing a flare in the light curve. In the framework of two- 
components jet the flare is produced by the deceleration of a 
moderately relativistic jet component as it interacts with the 
external medium. Patchy jets are characterized by large en- 
ergy fluctuations in the angular direction and can be imagined 
as multi-component jets, thus they can produce several bumps 
in GRB light curves. However, these bumps are typically too 
shallo w to explain the fa s t rise and decay observed in X-ray 
flares dZhang et al. ll2006h . lFan~& Wei I d2005h have shown that 
a FS-RS scenario can also explain the appearance of bumps 
in GRBs light curves; for appropriate model parameters, syn- 
chrotron RS emission can dominate in the X-ray band, pro- 
ducing a signature in the light curve. In the framework of ES, 
the interaction of the fireball with clouds of sufficiently small 
radii can produc e high variable GRB l ight curves and a delayed 
flaring ac t ivity dDermer et al. 2000; Dermer I l2007l) . Finally, 
Giannios (2006) showed that magnetic dissipation during ex- 
ternal shocks can produce fast-evolving and energetic flares if 
the dissipation regions have appropriate dimensions. 

Within group 2 models, X-ray flares c an be pr oduced by 



2005; Galli & Piro 200' 



LIS ( Burrows et al. 112005] : Wu et all 120 06) or ES JPiro et al. 



Panaitescu 1 120071) . The prolonged 



central engine activity could be due to a lon g-lived accretion 
of blobs of materi al on the central black hole (IKing et a. 1l2005l 



Pern a et al. 2006). In the context of LIS, X-ray flares are con- 
sidered prompt emission, that lasts for hundreds to thousands of 
seconds. We note that, in such a case, late internal shocks have 
to be tuned to produce X-ray emission rather than gamma-ray 
emission. In the case of ES by a long duration central engine 
the flare represents the onset of afterglow emission. This con- 
dition is verified if the central engine releases a significant frac- 
tion of the energy at late times, and when its duration is longer 
than the deceleration time (the so called thick shell case). In 
this case the onset of afterglow emission is delayed by several 
hundreds of seconds; we refer to this model as delayed external 
shock (DES). Recently. iPanaitescu I d2007l) has shown that X- 
ray flares could also be produced through the up-scattering of 
the forward shock emission by a more relativistic and lepton- 
enhanced shell emitted at the time of the flare appearance. 

X-ray flares appear mostly during the first phases of after- 
glow emission; consequently, independently on their origin, the 
delayed X-ray flare photons are expected to interact with the 
afterglow electrons by inverse Compton (IC) giving rise to de- 



layed high energy counterparts. In this work we study flares 
in the context of DES and assume an external medium with a 
smooth density profile (i.e a uniform interstellar medium (ISM) 
and wind-like medium). We decide to study flares in the context 
of ES, motivated by the spectral similarities observed in some 
bursts between X-ray flares and afterglow emission, which can 
be straightforwardly explained in this scenario. This scenario, 
depending on the spectral index of the electron population, can 
also account for those flares showing a spectral evolution of 
the order of 0.5-1.0, if the typical emission frequency is pass- 
ing through the observational band. Swift observations showed 
that several bursts present multi ple flares in their X-ray light 
curve, e.g. GRB 050730 dPandev et al.ll2006l) . In these cases, 
we can explain only one flare because only one flare can repre- 
sent the onset of afterglow emission, and the other flares have 
to be ascribed to the interaction of the fireball with a clumpy 
medium (in the framework of ES), or to any other of the mod- 
els of group 1 and/or 2. 

In the DES scenario the emission mechanism responsible 
for X-ray flares is synchrotron, and the flare photons can be 
IC up-scattered, thus producing flares in the GeV band. In this 
case, X-ray and high energy flares are produced by the same 
emitting region and electrons population, thus one expect that 
- depending on the energy - the two flares have similar tem- 
poral profiles and no significant delay. These predictions can 
be tested during the GLAST - Swift era and play an impor- 
tant role in order to discriminate between several models pro- 
posed to explain the origin of flares. Indeed, high e nergy flares 



can be produced also in the context of LIS mode ls (Wang et af 
20061 iFan & Piran Il2006t iFan et al. |[2007b . IWang et all d2006) 



have proposed two possible mechanisms that are able to pro- 
duce flares. In the first one, X-ray flares are produced through 
synchrotron emission in LIS, and the flare photons are self- 
IC scattered by the same electrons, producing the flare. In this 
case, the X-ray and the high energy flare come from the same 
region, thus they are expected to have a good temporal corre- 
lation as in the ES case. However, in the context of LIS, low 
Loren tz factors are required to produce long timescale X-ra y 



flares dFalcone et al. ll2005tlWang et al .1120061; IWu et al. 112006b 
and consequently less energetic flares are expected in the GeV 
band in comparison with the ES model. In the second mecha- 
nism, X-ray flares are produced by first-order IC emission from 
LIS, and these photons have a second order IC scatter with the 
afterglow electrons. X-ray flare photons need time to reach and 
scatter with the afterglow electrons. During this time the beam 
spreads out, and this causes a delayed and longe r high energy 
flare (for a detailed analysis see lFan et al. Id2007l) ). In addition, 
in this case we still expect high energy flares of lower intensity 
with respect to those expected by the ES. 

The paper is organized as follows; in Sect. [2] we study IC 
versus synchrotron emission in the so called standard case of 
a thin shell fireball. We first follow an analytical approach that 
can be applied during the fireball deceleration phase (Sect. 12. lV 
We evaluate the effect of pair production on the emitted radia- 
tion, and study the detectability of the IC component above the 
GeV band by the Large Area Telescope (LAT) aboard GLAST. 
We also develop a numerical model that permit the study of the 
fireball at times preceding the deceleration phase (Sect. l2~4i i. In 
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Sect. [3] we show that an ES by a thin shell fireball with syn- 
chrotron plus IC emission is able to model the delayed high 
energy emission observed in GRB 940217. We study the prop- 
erties of X-ray flares in the context of the DES (thick shell fire- 
ball), and determine the burst properties improving the detec- 
tion of high energy flares, see Sect. [4] In Sect. I4.2l we first pro- 
vide analytical formulas that allow one to predict, under rea- 
sonable assumptions, whether an X-ray flare with given ob- 
served properties can produce a detectable high energy flare. 
We then present a numerical study of X-ray and high energy 
flares by a DES (thick shell fireball), and apply this model to 
the individual cases of XRF 011030 and GRB 011121 (Sect. 
|5]l. We present our conclusions in Sect. [6] 



is the minimum post-shock electrons Lorentz factor, yi imax is 
the maximum electrons Lorentz factor and hv x is the peak en- 
ergy of the synchrotron em ission. Using the prescriptions of 
Panaitescu & Kumar I d2000h we find that the condition that IC 
emission occurs in the Thompson regime translates into: 



e eJim < 0.28 



f— )" 

\lkeVl 



(1) 



either for a fireball expanding in a uniform or a wind-like 
interstellar medium. The relative importance of IC and syn- 
chrotron emission is given by the Compton parameter Y, de- 
fined as the ratio of IC and synchrotro n luminosities L syn and 
L/c, see Eq. (3.1) of ISari & Esinl d200lb : 



2. Inverse Compton emission for a standard thin 
shell fireball in the framework of external shock. 

The emission mechanism explaining GRB prompt and after- 
glow radiation is synchrotron emission from relativistic elec- 
trons accelerated in internal and external shocks, respectively. 
We assume that the electrons are accelerated to a single rela- 
tivistic power law distribution of index p, i.e. N(y) oc y~P, with 
electron minimum energy y, = (m p /m e )e e (r - 1), where T is 
the fireball Lorentz factor and e e is the fraction of the fireball 



Y - Lic/L s 



(2) 



energ y going into relativistic electrons (Pan aitescu & Kumar 



2000). Photons produced by synchrotron emission can interact 
with the electrons by IC (or self synchrotron Compton: SSC), 
giving rise to a high-energy emission component. 

The synchrotron emission spectrum is described by four 
power law segments separated by the three characteristic fre- 
quencies, the absorption frequency v a , the injection frequency 
v, and the cooling frequency v c , during both the slow and fast 
cooling regimes. The IC component has a spectrum similar to 
that of the synchrotron component with the characteristic fre- 
quencies boosted by a factor ~ y 2 : v a jc - 2y 2 v„, vuc. — 2y 2 vj 
and v c jc = 2y 2 v c in the Thompson regime. I Sari & Esinl (1200 lh 
have shown that at frequencies greater than vyc the flux has to 
be increased by a logarithmic term with respect to the power 
law approximat ion. For simplicity, in this w ork we follow the 
prescriptions of iPanaitescu & Kumar (2000) and adopt for the 
IC spectrum a power law approximation, which is also a more 
conservative choice. 

2.1. A comparative study of Inverse Compton versus 
synchrotron emission 

In this section we study the relative importance of IC and 
synchrotron emission as a function of the parameters of the 
Fireball Model, i.e. the energy injected into the fireball £53, the 
external medium density n or the parameter A* in the case of a 
wind-like medium, and the efficiencies e e and eg- We initially 
study synchrotron and IC emission at time scales larger than 
the fireball deceleration time when analytical formulas can be 
used (for a detailed numerical study see Sect. |2.4l i. 

We also require that IC emission occurs in the Thompson 
regime, i.e that the energy of target photons yj(hv) is smaller 
than the electron rest energy m e c 2 . This condition translates 
into 7, < y^mox = Fm e c 2 /(hv x ) ( Wang et al. 2006I) . where y, 



We compute Y according to Eq. 35 of lPanaitescu & Kumar 



(2000), correcting the (e e /eg) ratio for the radiative efficiency 
T] (the efficiency 77 is the fraction of energy radiated via syn- 
chrotron and IC): 



Y=~ 



,2(3 - 5) e B 



17+I 



a/2 



1 



(3) 



with s=0(2) for an ISM (wind), 77 = 1 during the fast 
cooling regime and 77 = (v c /v,) (2_p ^ 2 during the slow cooling 
regime. Y is proportion al to (?7e e /eg) 1/2 i f J]e e /eB » 1 and to 
(j]£ e /e B ) ifr)e e /e B < 1 dSari &Esinll200lh . 

Clearly, the model parameters that mainly determine the 
relative importance of IC and synchrotron emission are e e and 
eg- We first study the importance of IC emission as a function 
of e e and eg in the case of a fireball expanding in an ISM, and 
set the remaining model parameters to £53 = 1.0 , n = 5.0, 
p = 2.5, and fix the redshift to z — 1. In Fig. [T] we plot Y as 
a function of (e^/e^) and at t = 500 s. With the parameters 
above and an initial fireball Lore ntz factor rp=100, the fireba ll 
deceleration time is td ec ~ 250 s (Panaitescu & Kumar 2000), 
i.e at t — 500 s the fireball is at the beginning of the deceleration 
phase. 

The red solid line in Fig.[T]corresponds to Y — 2, the purple 
solid line to Y — 5 and the blue solid line to Y — 10, while the 
black solid line is the line of IC and synchrotron equal lumi- 
nosity. 

The black region in Fig. [T]is forbidden because it implies 
£e > 1 . We have also checked that the Thompson regime condi- 
tion (that we express through Eq. [TJ is verified for any choice 
of parameters. The dashed line in Fig. [T] separates the region of 
parameters space where the fireball is in the fast cooling regime 
from that where the fireball is in the slow cooling regime, and 
corresponds to values of and (ee/es) that satisfy the follow- 
ing equation: 



€ B = 7.96 x 10- 2 (£ 53 «r 1/4 (l + Y)- l/l (e e /e B ) 



-1/2/ 



1/4 

(4) 

where T4 is the time expressed in unity of 1 day. This equa- 
tion tells us that for a given set of model parameters the size of 
the fast cooling region decreases with time, and that at a given 
time the higher the values of £53 and n are, the larger the size of 
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E M =1.0, n=5.0, p=2.5, z= 1 .0, t=500 s 




10 100 1000 



Fig. 1. Regions in the (ee/es), €b parameters space where the 
IC component dominates over the synchrotron component at 
t = 500 s for a fireball expanding in an ISM. The other model 
parameters are fixed to £53 = 1.0, n = 5.0, and p = 2.5. The 
redshift is z — 1. The black solid line corresponds to Y = 1, 
the red solid line corresponds to Y = 2, the purple solid line 
to Y = 5 and the blue solid line to Y = 10. The dashed line 
separates the region of the parameters space where the fireball 
is in the fast cooling regime from that where it is in the slow 
cooling regime. 

the fast cooling region is. For a given eg and {e e I^B) the relative 
importance of IC emission with respect to synchrotron emis- 
sion is greater during the fast cooling regime, thus IC emission 
is favored by high values of £53 and a medium with large den- 
sity n. The solid lines in Fig.Q]show that during the fast cooling 
regime, the relative importance of IC and synchrotron emission 
depends only on the (e e /e B ) ratio and increases with (e e /eB) for 
any values of e^. During the slow cooling regime 77 < 1.0, and 
to achieve the same value of Y greater values of (e e /e B ) are 
required. 

We also study how the vF v IC flux varies in the (e e /eB), 
£b parameters space (see Fig. [2j. We compute IC flux using 
Eq. |2] In this equation, synchrotron and IC luminosity should 
be calculated by integration over the total spectrum, but one 
can evaluate them and give order of magnitude estimates com- 
puting synchrotron and IC flux at the peak of the spectrum, 
i.e L syn =i v max F v (v max ) and L IC =2 v ICjnax F Vlc (v IC , max ), with 
v„iax = v« during the fast cooling regime and v max = v c during 
the slow cooling regime. For a fireball expanding in an ISM we 
find: 

(vF v ) IC = 2.4 x 10- 10 Z)^53en(— Vl + YT l Y 

W (5) 

T d '(1 + z) erg cm~ 2 s~ x 
during the fast cooling regime and 



E S5 = 1 .0, n = 5.0, p=2.5, z=1 .0, t=500 s 




10 100 1000 

e./e B 



Fig. 2. IC flux in the (tp/eg), eg parameters space at t — 500 s 
for a fireball expanding in a ISM. The other model parameters 
are the same as Fig. Q] The dashed line separates the region 
of the parameters space where the fireball is in the fast cooling 
regime from that where it is in the slow cooling regime. The red 
solid curve refers to an IC flux F IC — 5 x 10~ 10 erg cm~ 2 s~ x , 
the purple solid line refers to F IC = 5 x 10~ 9 erg cmr 2 s~ x , and 
the blue solid line refers to F IC = 10~ 8 erg cnT 2 s~ l . 

(vF v ) IC = 3 x l^Df 2i E^I\\ + Y)-^y4~T _ 

V £b I (6) 

Tf\\+zf A erg en 2 s- 1 

during the slow cooling regime for p = 2.5. During the fast 
cooling regime the IC flux increases with e e . It does not de- 
pend on either e B nor on the density n of the external medium. 
This implies that the estimates presented in Fig. [2] are not 
strongly affected by the uncertainties regarding the external 
medium density n. During the slow cooling regime the ra- 
diative efficiency r\ is proportional to n( p-2 " 2 (l + Y) (p - 2) and 
Y ~ n (p-2)/2(4-p) > thus f rom Eq.|6]we find (vF v ) IC ~ n { P- 2)l{A -P\ 
i.e IC flux increases with n. 

In the case of a fireball expanding in a wind-like medium 
we find similar results, i.e the Compton parameter Y and the 
IC peak flux have similar behaviors as a function of the model 
parameters and similar values to those found in the ISM case. 
The fast-to-slow cooling transition occurs for: 

e B = 6.6 x 10- V 2 (l + F)- 1/2 (^) _1/ V] /2 (l + zT 1 ' 2 (7) 

This equation tells us that for a fixed time, the higher the 
wind density parameter A, is, the larger the region of the pa- 
rameters space where the fireball is in the fast cooling regime is, 
i.e. IC emission is also favored by large densities when the fire- 
ball expands in a medium with a wind density profile. Contrary 
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to the ISM case, the extension of the fast cooling region does 
not depend on the fireball energy £53. Using the same approach 
followed for the ISM case, we find that for a fireball expanding 
in a wind medium density profile the IC flux is: 



(vF v ) IC = 4.7 x 10- 10 D^£ 53 e B ( J)(l + Yy x Y 



T d l {\ + z) erg cm 2 s 1 



(8) 



T yy — 1 the pair attenuation cannot be neglected, and the spec- 
trum is limited by a cutoff energy E cut : 



0.5E 



-1/4 -ygi 

n 4 '' 



2tg-l) _ (p-2) 



e e-i £ b-2 V+Y) Vp T d 4p (1+z)— TeV 

(13) 

Similarly, for a fireball expanding in a wind-like medium 
we find the cutoff energy to be: 



during the fast cooling regime and 



<yF v ) IC = 7.2 x 10- 9 D u 2 8 £ 53 Ay 2 4g) 3/2 (l + Y)^ 2 YT^ 2 



(1 + z) 3/2 erg cm 2 s 1 



(9) 



during the slow cooling (for p = 2.5). For a fireball ex- 
panding in a wind-like medium during the fast cooling regime, 
(vF v )ic increases with e e , while it does not depend on eg and 
the wind parameter A«, as for the ISM case. During the slow 
cooling regime the radiative efficiency 77 is proportional to 

thus from Eq. [9] we 



A!f" 2) (l + Y) ( p- 2) and Y ~ A 



(p-2)l(A-p) 



find(vF v ) /c ~A^- 2)/(4 - rt , 
(vF v )/c increases with A,. 



i.e. during the slow cooling regime 



2.2. Pair attenuation estimates 

Here we estimate the attenuation of IC emission by pair pro- 
duction due to yy interactions. This process could be very im- 
portant because it causes internal absorption which determines 
a cutoff in the spectrum of the emitted radiation at high ener- 
gies. Given a photon of energy E y , the threshold energy E,i, for 
the production of an electron-positron pair is: 



2(m f c 2 ) 2 r 2 
E y (l - cosB) 



(10) 



where 6 is the impact angle between the two photons 
in the shell frame. T a king i nto account the prescriptions of 
Panaitescu & Kumar ( 2000h . for a fireball expanding in an 
ISM we find that the threshold frequency v f /, where pair pro- 
duction becomes important is: 



v, h = 2.5 x 10 13 £ 5 ^ n 



i5 F i/4 M -i/4 T -3/4 



(1+z) 



-5/4/ E y \ 
\lTeV) 



Hz (11) 



The optical dept h to pair productio n for a photon of energy 



E y is (see Eq. (2) of Wang et al. (2004): 



/ 11 \ o-TDl m F v (v th ) 
1 180/ 4T 4 c 2 hpt 



(12) 



where (11/180 )ct7- is an analytical approximation for the 
yy cross section ( ISvensson 1 119871) . D[ um is the source dis- 
tance and p is the spectral index of the radiation emitted by 
the source. During the first phases of afterglow emission, i.e. 
around ~ 500s, for typical GRB model parameters v c < v, < v,;, 
thus F y (v th ) = F Vtmax (v th /v i )- p/2 (v c /v i ) 1/2 and/3 = p/2. When 



0.9F^ P A~ 



2(p-l) lp-2) 

" - lp (1+Y) 2/P T, 



' (1+z)" 



TeV 
(14) 

For typical model parameters E cut ~TeV, around 500 s, both 
for a fireball expanding in an ISM and in a wind-like medium; 
consequently our model predictions (which do not account for 
the pair attenuation) should not in general be affected by in- 
ternal absorption up to energies of the order of TeV. We note 
also that for typical model parameters around 1 day the cutoff 
energy is still above the LAT energy band. 

2.3. Detectability of the Inverse Compton component 
in the GeV band with the GLAST Large Area 
Telescope 

Now we compare the IC flux predicted in the context of the ES 
by a thin shell fireball with the GLAST Large Area Telescope 
(LAT) sensitivity during both early and late afterglow emission. 
We estimate the LAT sensitivity adopting the same criterion of 
Zhang & Meszaros I (120041) . that a 5<x detection is made when 
at least 5 photons are collected if the instrument is source dom- 
inated. We evaluate the LAT sensitivity between 100 MeV and 
200 GeV, the spectral range where the instrument point spread 
function (PSF) and effective area A e ff are well characterized. 
We integrate at each energy E the burst spectrum in an energy 
range AE = E/2, and also take into account the time when 
the LAT becomes background dominated. When the instrument 
is source dominated the LAT sensitivity F,/, decreases linearly 
with time according to the following law: 



F th (E) = 



A eff (E)t 



cm 2 s 1 



(15) 



where E is the photon energy and t is the inte- 
gration time. The most important contribution to the 
background comes from the observed high latitude dif- 
fuse flux, ~ 1.5 x 10 _5 cm _2 s _1 sr _1 above 100 MeV, see 
http://www-glast.slac.stanford.edu/software/IS/glastJat_performance.htm 
When the instrument is no longer source dominated its sen- 
sitivity starts to decrease with f 1 ^ 2 . The sensitivity curve is 
displayed with a black dot-dashed line in Figs. QT| and Q~3] 
for a short integration time (500 s), and in Figs. [5] and [7] for 
a long integration time (10000 s). We also evaluate the LAT 
sensitivity threshold as a function of the integration time at 1 
GeV, the observational energy where we display the predicted 
light curves, see Sects. [2~4l and |4. 21 At 1 GeV the on-axis LAT 
effective area is A e ff 0.85 m 2 and the PSF is of the order of 
1°, thus the LAT sensitivity is ~ 1.9 x 10~ 9 erg cm~ 2 s~ x for an 
integration time t = 500 s, and is ~ 9.4 x 10~" erg cm~ 2 s~ l 
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for an integration time t = 10000 s, and the LAT is source 
dominated up to an integration time of ~ 61500 s. 

We compute the IC flux for a burst located at redshift z=l 
at the ob servational energy of 1 GeV (following the prescrip- 
tions of Panaitescu & Kumar (EoOOh . In Fig. [3] we compare 
our predictions with the 1 GeV LAT sensitivity (grey solid line) 
for an integration time of 500 s. It shows that LAT will be able 
to detect IC emission only during the fast cooling regime, and 
for a very narrow region of the (e^/eg), €b parameters space. 
However, if the integration time is sufficiently long (but the in- 
strument is still source dominated) LAT will be able to detect 
IC emission for a significantly larger region of the parameters 
space, see Fig. [4] We find similar results for the detectability 
of the IC component associated with a fireball expanding in a 
medium with a wind density profile; also in this case the de- 
tectability of the IC component increases with the integration 
time if the LAT is source dominated. 



0.0100 F 



0.0100 



<S 0.0010 



0.0001 
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1000 



Fig. 3. vF v of the IC at 1 GeV in the (^/eg), £b parameters 
space at t — 500 s for a fireball expanding in an ISM. The other 
model parameters are the same as in Fig. Q] The dashed line 
separates the region of the parameters space where the fireball 
is in the fast cooling regime from that where it is in the slow 
cooling regime. The grey solid line represents the LAT sensitiv- 
ity for an integration time of 500 s. The purple solid line refers 
to an IC flux (vF v ) IC = 5 x 10~ 10 erg cm~ 2 s~ x , and the red solid 
line to (vF v ), c 



2x 10~ 10 erg cm- 2 s' 1 . 



We can make similar considerations for the detectability of 
the IC component by the gamma-ray imaging detector (GRID) 
aboard AGILE (Astro-rivelatore Gamma a Immagini LEggero; 
launched April 2 3, 2007), which is s ensitive in the range 30 
MeV-50 GeV (see Tavani et al. ( 2006h for a recent review). To 
make a comparison between the GRID and the LAT sensitivity 
it is important to know the time when the GRID becomes back- 
ground dominated. At early times the GRID is source domi- 



0.0010 




SLOW COOLING 



0.0001 



10 



100 
e./e B 



1000 



Fig. 4. vF v of the IC at 1 GeV in the (e e /es), 6b parameters 
space at t — 10000 s for a fireball expanding in an ISM. The 
other model parameters are the same as in Fig. [T] The dashed 
line separates the region of the parameters space where the fire- 
ball is in the fast cooling regime from that where it is in the 
slow cooling regime. The grey solid line represents the LAT 
sensitivity for an integration time of 10000 s. The red solid line 
refers to an IC flux (vF v )jc = 2 x 10~ 10 erg cm~ 2 s~ l , and the 
purple solid line to (vF v )i C = 5 x 10~ 10 erg cm~ 2 s~ l . 



nated and we can evaluate its sensitivity using Eq. Q3] as for 
the LAT. At late times the GRID is background dominated and 
its sensitivity decrea ses as r 1 ^ 2 . We take the GRID sensitiv- 

d2006l) and we find that it 



ity at 50 hours from iTavani et al 



lies above the extrapolation of eq. Q3] This implies that at 50 
hours the GRID is background dominated. The transition from 
one regime to the other takes place around 2000 s. The GRID 
on axis effective area at 1 GeV is ~ 550 cmr 2 which implies 
that its sensitivity is ~ 15 times lower than that of the LAT, 
i.e ~ 2.9 x 10~ 8 erg cm~ 2 s~ x for an integration time of 500 s. 
For longer integration times, when the GRID becomes back- 
ground dominated, its sensitivity goes as F,k = 5.5 x 10 _7 i -1 ' 2 , 
which gives a flux threshold of ~ 5.5 x 10~ 9 erg cm~ 2 s~ l for 
an integration time of 10000 s. A comparison of GRID sensi- 
tivity with Figs.[3]and|4]shows that it will not be able to detect 
ES IC emission for a burst located at z=l for both short and 
longer integration times. This is consistent with the findings of 
Gou & Meszaros I d2007l) . 



2.4. Inverse Compton spectra and light curves 

Here we would like to asses in detail the maximum contribution 
of IC emission and how it varies with time in the framework 
of standard ES (i.e thin shell fireballs). Analytical equations 
are valid only at times larger then the fireball deceleration time 
and permit us to model only the fading portion of the burst light 



A. Galli & L. Piro: GRB high energy emission 



7 



curve. Now we adopt a numerical approach whi ch allow us to 
describe also the rising part of the light curve (see lGaili & Piro 
d2006l) for details). 

We start with the case of a thin shell fireball expanding in 
an ISM, which remains in the fast cooling regime during a large 
time interval of afterglow emission, e.g. up to thousands of sec- 
onds after the burst, because this is the condition that maxi- 
mizes the IC component. In Sect. 12. II we have shown that dur- 
ing the fast cooling regime the relative importance of IC and 
synchrotron emission is determined only by the (e e /eB) ratio, 
and the IC flux increases linearly with e e while substantially 
does not depend on e B . We thus assume a large (but reasonable) 
value of e e that also satisfies the request that IC emission occurs 
in the Thompson regime (see Eq. [TJ, and vary the value of e B . 
We choose the other model parameters consistently with the re- 
quests of a thin shell and fast cooling regime. Equation |4] tells 
us that the time the fireball remains in the fast cooling regime 
increases with £53 and n, and Eq. tells us that IC emission 
increases with £53; this motivates us to assume large values 
of fireball energy and external medium density such as £53 = 1 
and n — 5. In an ISM the thin shell condition reads (IPiro et al 



2005): 



1/3 



r 8/3 < 6 2 

1 , ^ u.z. 



(16) 



where T2 is the prompt emission duration in unity of 100 s 
and ro,2 is the initial fireball Lorentz factor in unity of 100. We 
fix To = 120, which ensures the thin shell condition is satisfied 
for a prompt emission duration as long as ~ 200 s. We note 
that the value of To affects only the first phases of ES emission, 
while during the deceleration phase, neither synchrotron nor 
IC emission depend on To. Finally, we fix the redshift z and the 
spectral ele ctron population index p to canonical values as z=l 
and p-2.5 (IDe Pasquale et al. 1 2006). 

We present the effect of (e e /e B ) on synchrotron and IC 
emission in Fig. [5] where we give spectra at ~ 10000 s after 
the burst. We evaluate the IC spectrum cutoff energy due to 77 
interactions using Eq.[13]for each set of model parameters, and 
display it with a thin vertical solid line. We find that in this 
case, pair attenuation has little to not effect on the predicted IC 
emission in the LAT energy band. The black dot-dashed line 
represents the LAT sensitivity obtained by applying the method 
explained in Sect. 12. 3 1 for an integration time of ~10000 s. We 
notice that with standard model parameters, LAT can detect 
afterglow IC emission for a sufficiently long integration time 
(see also Fig.|6]l. These spectra also show that as expected, the 
smaller the value of s B is, which corresponds to greater val- 
ues of the (ee/eg) ratio, the higher the IC emission with respect 
to the synchrotron emission is (see Sect. |2). At early times IC 
emission dominates over synchrotron emission only at high en- 
ergies, i.e., in the GeV band, but at late times IC emission dom- 
inates on synchrotron emissi on also at lower ener gies, i.e. in the 
hard X-ray/soft-7-ray band (ICorsi & Piro 112006b . In Fig. [6] we 
show how synchrotron (dashed line) and IC (solid line) light 
curves vary with e B and compare these light curves with the 
LAT threshold (black dot-dashed line). 
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Fig. 5. Synchrotron (dashed lines) and IC (solid lines) spectra 
for a thin shell fireball expanding in an ISM at t ~ 10000 s. 
The model parameters are Fo=120, £53 = 1.0, n = 5 and e e 
0.2. Green, blue and red spectra are obtained respectively for 
e B = 0.0001, e B = 0.001, and e B = 0.01. The thin solid lines 
indicate the cutoff energy due to pair production for the three 
IC spectra. The black dot-dashed line is the 5 cr LAT sensitivity 
for an integration time of ~ 10000 s. 
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Fig. 6. Synchrotron (dashed line) and IC (solid line) light 
curves for a thin shell fireball expanding in an ISM. The model 
parameters are the same as in Fig. [5] Synchrotron light curves 
are computed at 10 18 Hz and IC light curves are computed at 
2.4 x 10 23 Hz. The green curves correspond to e B = 0.0001, the 
blue curves correspond to e B = 0.001, and the red curves cor- 
respond to £ B = 0.01. For these choices of model parameters 
at ~ 500 s the fireball is in the fast cooling regime. The black 
dot-dashed line shows the 5 cr LAT threshold at 1 GeV as a 
function of the integration time. 



If we change the model parameters so that during the 
first hundred seconds of ES emission the fireball is still in 
the slow cooling regime, in this case the 1 GeV IC flux is 
~ 2x 10~ 12 erg cmr 2 s~ l at around 10000 s after the burst, while 
the LAT sensitivity at this time is ~ 9.4 x 10~ n erg cmT 2 s~ x . 
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Thus, in such a case, there is no possibility for the LAT to detect 
the IC component. 

Now we discuss our predictions for a fireball expanding in 
a wind-like medium. In this case we require that the fireball 
remains in the fast cooling regime up to thousands of seconds 
after the burst, to maximize IC emission. We chose the model 
parameters according to the same criterion followed in the ISM 
case. During the fast cooling regime Y depends only on the 
(e e /efl) ratio and the IC flux increases with the fireball kinetic 
energy £53 (see Eq. [8) as for the ISM case. Equation [7j shows 
that the time the fireball remains in the fast cooling regime in- 
creases with the wind parameter A*, and Eq. [8] indicates that 
IC emission increases with £53. In order to have a significant 
high energy emission we thus assume £53=1 and A, = 0.5. In 
a wind the thin shell condition reads dPiro et al. I l2005l) : 



T 2 A*ri 



0.2 



£53 



<0.06 



(17) 



In this case we adopt To = 50, so that with the parameters 
above the thin shell condition is verified for a prompt emission 
duration < 200 s. 

In Fig. |7] we show synchrotron and IC spectra as a function 
of eg at t ~ 10000 s: as expected, the relative importance of IC 
and synchrotron emission increases with decreasing eg value, 
i.e with higher values of (e e /es). Most of all this figure shows 
that, as for the ISM case, in this case LAT can detect afterglow 
IC emission if the integration time is sufficiently long (see also 
Fig.©. 
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Fig. 7. Synchrotron (dashed lines) and IC (solid lines) spectra 
for a thin shell fireball expanding in a wind at t ~ 10000 s. 
The model parameters are To=50, £53=1, A*=0.5, and e e -0.2. 
Green curves correspond to eg - 0.0001, blue curves corre- 
spond to £b = 0.001, and red curves correspond toes = 0.01. 
The thin solid lines indicate the cutoff energy due to pair pro- 
duction for the three IC spectra. The black dot-dashed line is 
the 5 <x LAT sensitivity for an integration time of 10000 s. 

When the fireball is already in the slow cooling regime 
during the first phases of fireball evolution, as happens for a 




Fig. 8. Synchrotron (dashed line) and IC (solid line) light 
curves for a thin shell fireball expanding in a wind. The model 
parameters are the same as in Fig. [7] Synchrotron light curves 
are computed at 10 18 Hz and IC light curves are computed at 
2.4 x 10 23 Hz. The green curves correspond to ££=0.0001, the 
blue curves correspond to eg =0.001 and the red curves cor- 
respond to t?B=0.01. The black dot-dashed line shows the 5<t 
LAT threshold at 1 GeV as a function of the integration time. 



fireball expanding in an ISM LAT will again be unable to de- 
tect the IC component, even for sufficiently long integration 
times. In fact, the predicted 1 GeV IC flux at ~ 10000 s is 
~ 3 x 10~ 13 ergcrrT 2 s~ l , i.e. well below the corresponding LAT 
sensitivity. 

In Sect. l2.3l we have shown that the GRID sensitivity is ~ 15 
times lower than that of the LAT when both the instruments are 
source dominated, and also lower for longer integration times 
(such as 10000 s) when the GRID is background dominated 
while the LAT not yet. As expected (see Sect. 12.31 1 we find that 
at any time the predicted flux is below the AGILE sensitivity. 



3. GRB 940217 in the framework of the external 
shock model 

GRB 940217 was an energetic and long burst observed by the 
Energetic Gamma-Ray Experiment Telescope (EGRET) above 
30 MeV with a duration of ~ 5400 s a nd a possible emis - 
sion of an 18 GeV photon around 4500 s dHurley et al. II 1994 . 
The origin of this delayed high energy emission is still de- 
bated. At such high energies, IC is one of the best candi- 
date emission mechanisms. We investigate whether IC emis- 
sion from the afterglow of a "standard " thin shell fireball 
could expl ain the delayed high energy emission occurring in 



this event. Hurley et al. I (119941) showed that above 30 MeV 
the best fit of the mean spectrum of the delayed high energy 
emission is a power law with photon index y = 2.83 + 064 
(the 18 GeV photon is consistent with this fit at the 99% 
confidence level). The integration of this power law gives a 
fluence S = 7 x 10~ 6 erg cm" 2 , which implies a mean flux 
F 



1.4 x 10 9 erg cm 2 s l . Figure [2] shows that a flux as high 
as that observed for GRB 940217, and greater also, is achieved 
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for typical GRB model parameters in the case of thin shell fire- 
balls when the fireball is in the fast cooling regime. 

During the fast cooling regime the photon index of the IC 
spectrum is y= 3/2 for v c jc < v t s < Vijc, and y = (| + 1) 
for < v obs- If) at the ti me of the delayed emission, 

1S below the observational band, for p-2.5 one obtains 
y=2.25, which is c onsistent with the best fit value found by 



Hurley et al. 



(119941) . We thus require that at the time of the 
delayed emission vijc < 30 MeV. As noted by IWei & Fan 



(2007), the observational data indicate that the flux of the de- 
layed emission is roughly constant with time. We succeed in 
finding a solution that accounts for both the spectral and tem- 
poral properties of the delayed high energy emission and its 
mean flux, as presented in Fig.s|9]and[T0] In particular in Fig. 
TU] we compare the observed spectrum (green points) with that 
predicted in the framework of ES at ~ 500 s (red-dashed line) 
and ~ 5000 s (black solid line). We note that around 5000 s the 
predicted spectrum is fully consistent with the requirement that 
Vcjc < v ijc <30 MeV. Around 500 s the injection frequency 
Vjjc is still above 30 MeV, however the predicted spectrum is 
still in satisfactory agreement with the observed one. Below 
v c jc the IC flux goes as r' 3 and for v c jc < v < Vijc it goes as 
f 1 ^ 8 . This explains why the predicted flux is roughly constant 
along the delayed high energy emission timescale. However, 
we note that such a model is unable to account for the delayed 
18 GeV photon. 



10-' 



Time [sec] 

Fig. 9. IC light curve of GRB 940217 in the 30 MeV-500 MeV 
energy range in the context of ES for a thin shell fireball ex- 
panding in an ISM. The model parameters are £53 = 5.0, 
n = 3.0, e e = 0.07,e B = 0.001, p = 2.5, and the redshift is 
Z = 1 . IC flux is roughly constant up to thousands of s, consis- 
tent with the observed delayed high energy emission. 




Fig. 10. IC plus synchrotron spectrum of GRB 940217 at 500 
s (red dashed line) and 5000 s (dark solid line) after the burst 
in the context of ES. The model parameters are the s ame as in 
Fig-El Data points are taken from lHurlev et al. I (11994). 



4. Prediction of GeV flares in the framework of 
DES 

4.1. X-ray flare properties 



We se lect from the sample of "Gold" flares of iFalcone et al 
(2007) (i.e. the sample of X-ray flares with better statistics) a 
sub-sample of X-ray flares appearing in the first ~ 1000 s af- 
ter the burst and compare their spectral properties with those of 
the underling X-ray light curve. Flares and light curves are both 
fitted with a power law model. We retrieve the spectral p roper- 
ties of the X-ray continuum from lWillingale et al. I d2007l) . who 
have performed the temporal and spectral analysis of a large 
sample of Swift GRBs (~ 100) in terms of two components, the 
first one attributed by them to the prompt y-ray emission and 
the initial X-ray decay, and the second one thought to be due to 
the external shock, and developing in the afterglow. We list our 
sub-sample of flares in TableQ] where we give th e flare number, 
the en d-time of the flare t en d (as determined by Falcone et al. I 
d2007h ). the power law photon index of the flare Tfi are , that of 
the prompt decay F p d and that of the afterglow decay T a d, the 
time T a at which the power law afterglow decay starts, and the 
difference Ar = Tfi are -Y a d between the flare and the afterglow 
decay power law photon indices; all the quantities reported in 
this Table are at 90% confidence level. We note that T a is al- 
ways larger than the time of the flare appearance. We find that 
for 16 flares out of 28, the flare photon index is not consistent 
with that of the prompt decay, nor with that of the afterglow de- 
cay. For 5 out of 28 flares Tfi are is consistent with both T p d and 
F fl rf. For one flare Y fi are is consistent only with T pi j, and for 5 
out of 28 flares Tfi are is consistent only with r a( j. In total, there 
are 10 out of 28 flares for which the spectrum is consistent with 
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that of the second component of GRBs light curve (that can be 
ascribed to the external shock). We also note that in the ma- 
jority of cases the spectral difference between the flare and the 
afterglow decay spectrum is AF ~0.5. Such hard-to-soft spec- 
tral evolution can be easily accounted for in the context of the 
ES model. Indeed, during the fast cooling regime, when the in- 
jection frequency v m passes into the observational band, there 
is a hard-to-soft spectral evolution of Ar = (p — l)/2, with 
0.5 < Ar < 1 for 2.0 < p < 3.0. Similarly, during the slow 
cooling regime, when the cooling frequency v c passes into the 
observational band there is an hard-to-soft spectral evolution 
of AF =0.5. If the fireball expands in a wind-like medium, the 
cooling frequency increases with time, thus one can also ex- 
plain a soft-to-hard spectral evolution of Ar =0.5, that was ob- 
served in some cases (see TableQ]). 

Another important quantity is the fireball kinetic energy £53 
that significantly affects the predicted synchrotron and IC flare 
emission (see e.g. Eqs.|5]and[8]); this quantity can be estimated 
from the emitted isotropic energy Ei S0 of the burst correcting 
it for the fireball radiative efficiency rj. In Sect. 12.41 we have 
shown that a favorable condition for a detectable IC high en- 
ergy emission is to have a burst with £53 >1. It is thus im- 
portant to understand what fraction of GRBs satisfy this en- 
ergy requirement. T o this end we have further selected from 



Falc one et al. J2007h a sub-sample of Swfit GRBs with known 



redshift, and report them in Table [2] Wh en it i s possible we 
take the value of E-, so directly from lAmati I (12007). For the other 
bursts we estimate Ej SO according to the following method. The 
burst isotropic energy E, so is given by the integra l of the burst 
-10 4 keV band (see 



Amati et al. 



spectr al model in the 1 

(Eoolh . In this energy band the spectrum of a burst is typically 
described by a Band model ( Band II 1993b . thus in order to cal- 
culate this integral one needs to know the two power law pho- 
ton indeces and the peak energy of the Band function. The only 
information we can retrieve from the on-line Swift GRBs table 
Q are the 15-150 keV fluence S and the prompt emission pho- 
ton index F of the burst measured by the Burst Alert Telescope 
(BAT). We therefore run a set of integrations by varying E pen k 
(for given values of spectral indices, see below) and derive 
the corresponding _?; so using the 15-150 keV BAT fluence to 
normalize the spectrum. The chosen value s of E oea k and Ej S p 
are those sati sfying the Amati relationship dAmati et al. 1120021 



AmatD|200j_b. In each integration, depending on the value of 



Epeak, we identify the observed power law photon index with 
one of the two indices of the Band function, fixing the other 
one to a canonical value (-1 and -2 for the low and the high en- 
ergy power laws respectively). We find that ~ 30 % of the bursts 
of our sample have isotropic energy E, so > 10 53 erg (see Table 
12). Finally we assume a radiative efficiency 77 ~ 0.1, which im- 
plies that ~ 30 % of the GRBs with flares in their X-ray light 
curve have fireball kinetic energy > 10 54 erg. 



4.2. Inverse Compton and synchrotron emission in the 
delayed external shock scenario: X-ray and high 
energy flares 

In the preceding section we have summarized some of the more 
relevant flare properties and have shown that there are flares 
for which temporal-spectral properties can be potentially ex- 
plained in the framework of external shocks, and more specif- 
ically, in the framework of th e DES scenario (i.e. thick shell 
fireballs, ( I Sari & Piran Il999h l. In particular, we have shown 
that a DES can also potentially explain spectral variations of 
the order of 0.5-1.0, depending on the radiative regime and 
electron population spectral index. In the DES scenario the du- 
ration of the central engine activity t eng is larger than the time 
tdec when the fireball starts to decelerate, and most of the en- 
ergy is transferred to the surrounding material at late times, 
around the end of the engine activity. The flare would thus be 
produced by an ES caused by a central engine whic h remains 



active until the time of the flare occurrence (see iPiro et al 



d2005l) : iGalli & Piro I d2006l) for details). We also stress that in 
this context, with the flare being ascribed to the onset of af- 
terglow emission, only one flare can be explained. To explain 
multiple X-ray flares other mechanisms, such us LIS by a long 
duration central engine, refreshed shock, or ES with a clumpy 
medium, are required. An important implication of models pro- 
posed to explain the origin of flares is that X-ray flare photons 
can be IC scattered, thus producing high energy (MeV to TeV) 
flares. In particular, in the context of DES the X-ray flare pho- 
tons are IC scattered by the afterglow electrons, and X-ray and 
high energy flares are produced by the same source region and 
electron population. 

Before studying in detail the properties of the predicted IC 
emission associated with flares, we first provide a rough esti- 
mation of the mean flux and the peak energy of the high energy 
flare predicted by the DES scenario as a function of X-ray flare 
flux and peak energy. This allows a first-order prediction of the 
detectability of high energy flares based solely on their lower 
energy properties. The IC flux can be estimated from that ob- 
served in X-ray th rough the following general equation (see 
Sari&EsinlfeOOhY ): 



VpjcFvjc = Yv p F v (18) 
where v p jc and v p are the peak frequencies of the IC and 
synchrotron components with F v jc and F v the related peak flux 
densities, and Y is given by Eq. [2] The equation above tells us 
that the flux, and consequently also the fluence, associated with 
the IC component is roughly Y times those of the X-ray com- 
ponent. |Fakojie_eialJ (12007) have derived the mean (0.2-10.0) 



http://swift.gsfc.nasa.gov/docs/swift/archive/grb_table.html/ 



keV flare fluence to be S sj3x 10~ 7 erg cmr 2 and the typical 
flare duration to be few hundred seconds for flares appearing 
within ~ 1000 s after the burst. This implies a mean X-ray flare 
flux of as 3 x 10~ 9 erg cm~ 2 s~ l , and a mean 1 keV X-ray flux 
density of « 0.5 mJy for a typical flare spectrum. On the other 
hand, the Compton factor Y strongly depends on the values of 
the parameters e e and eg. For example, at equipartition Y is of 
the order of unity, and in this case we thus expect that the flu- 
ence of the high energy flare is of the same order of magnitude 
as that of the X-ray flare, i.e. « 3 x 10~ 7 ergcmT 1 . Suchafluence 
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Table 1. Properties of power law spectral fits to the Gold sample of flares of Falcone et al. (2007) and to the total X-ray light 
curve. The spectral properties of the global X-ray light curve are taken from Willingale et al. " d2007h . Errors are at 90% confidence 
level. 
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Table 2. Sample of Swfit GRBs with known redshift having flares in their light curve. We report in this table their prompt 
emission prop erties, and in particular their isotropic energy Et so . When it is possible we take the value of Ej SO directly from 
Amati I (120061) : we mark these bursts with *. For the other bursts we assume the validity of the Amati relation, and determine the 
values of the intrinsic burst peak energy E pea k and Ej SO implied by the observed fluence of the burst itself. 
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is below the 5<x 1 GeV LAT fluence threshold, ~ 10~ 6 erg cmr 2 , 
and is of the same order of magnitude as the 5cr 100 MeV LAT 
fluence threshold, ~ 2.7 ~ 10~ 7 erg cmr 2 (with the LAT on- 
axis 100 MeV effective area ~ 3000 cm 2 ; see Sect. for the 
estimation of the LAT sensitivity). When Y is above unity (this 



is common during the fast cooling regime for typical model 
parameters; see Fig. [U, the fluence carried by the high energy 
flare can be larger than that carried by the X-ray flare, and the 
detection of high energy flares improves. It is important to note 
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that the detection of high energy flares strongly depends on the 
peak energy with respect to the LAT energy band. 

In Sect. 12.11 we have shown that we expect similar fluxes 
for a fireball expanding in an ISM and in a wind-like medium. 
Therefore, in this section we focus our attention only on the 
case of a thick shell fireball expanding in an ISM. We now esti- 
mate the peak energy of the IC component as a function of the 
observed X-ray flare flux and peak energy. To this end we first 
determine the values of the model parameters that satisfy the 
observed flare properties in the context of the DES scenario. 
We reasonably assume that at the time of the flare appearance 
the fireball is in the fast cooling regime, and that the injection 
frequency v m (i.e. the flare peak frequency in the v - vF v space) 
is below the observational frequency v b s - Under these assump- 
tions, if we fix the fireball kinetic energy to £53 = 10 (we have 
shown in the preceding section that this is a reasonable choice), 
p = 2.5 and z - 1.0, the flare mean X-ray flux and peak fre- 
quency v m depend only on the two parameters e e and eg, and 
we can search for a pair of e e and eg values accounting for 
the observed flare properties. When v t, s > v m , the synchrotron 
flux density is given by F v = F Vc (v hs/v m y p/2 (v c /v m y /2 . In this 
equati on we substitute F v „, v „ , and v r using Eqs. (64), (22) and 



(27) o flPanaitescu & Kumar (2000), respectively. We thus find 
two equations relating e e and eg to the 1 keV X-ray flare flux 
density F\k e v and peak frequency v m , and the time t b S when 
we measure the flare properties: 
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where Z?; SOj 53 is the emitted isotropic energy in unity of 10 
erg, and 771 is the fireball radiative efficiency in unity of 0.1. 
The ratio Ei S0 ^lrj is the way we estimate the fireball kinetic 
energy £53. At this point we can finally write the Compton fac- 
tor Y, the peak energy of the high energy flare v p jc an d its flux 
VpjcF v jc, in terms of only F\ u e v, and the flare peak energy v m 
and peak flux density F Ym : 
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The analytical equations given above hold only during the 
fireball deceleration phase. We thus take f o f, s =200 s in our cal- 
culations, i.e. t i, s is larger than the fireball deceleration time 



(which corresponds to the flare maximum in the DES scenario), 
but sufficiently close the peak time of the flare. From the pre- 
ceding Eqs. we find that if we observe an X-ray flare with 1 
keV flux density Fut e v=0.5 mJy and peak frequency v,„=10 eV, 
and the density of the uniform interstellar medium is n=10, for 
t b s ~ 200 s then Y ~1, We also expect a high energy flare with 
peak flux ~ 3 x 10~ 10 erg cmr 2 s~ l and peak frequency v p jc ~ 
25 MeV. If the X-ray flare peak frequency is v m =l keV and 
the other quantities do not vary, then Y ~100 and we expect a 
high energy flare with peak flux ~ 10~ 7 erg cm~ 2 s and peak 
frequency V p jc ~ 5 TeV. 

We now present two numerical solutions corresponding to 
two flares with 1 keV flux density of ~ 0.5 mJy and peak en- 
ergies around 10 eV and 1 keV, respectively. We take the den- 
sity of the external medium n — 10, the initial fireball Lorentz 
factor To = 150, and values of e e and eg that satisfy the con- 
ditions of fast cooling regime, Thompson regime and thick 
shell fireball (which can be obtained reversing Eq.fToTl. We as- 
sume a typical value of 500 s for the time of the fl are occur- 
rence dChincarini et al. 2007 : Falcone et al. 2007 ) and apply 
the thick shell c ondition, shifting the origi n of time to this in- 
stant dPiro et al. ll2005llGalli & Piro Il2006h . f = 500 s. 

We present in Figs. [TT] and [12] a numerical solution cor- 
responding to a flare peak energy of ~ 10 eV. Figure QT] shows 
synchrotron (solid line) and IC (dotted line) spectra around 200 
s after the start of the flare; the predicted high energy flare can 
be detected by the LAT up to ~ 1 GeV. In Fig.[T2]we give syn- 
chrotron (1 keV; solid line) and IC (100 MeV; dotted line) light 
curves, and compare the latter with the 100 MeV LAT sensi- 
tivity (dot-dot-dot-dashed line). When the peak energy of the 
X-ray flare is around 10 eV the peak energy of the high energy 
flare is around 200 MeV and there is a good temporal correla- 
tion between the X-ray and high energy flares; this is expected 
in the context of the DES scenario, when the peak energies of 
X-ray and high energy flares are below or very close to the rel- 
ative observational bands. 

Fig.sQ~3]and[T4]present a numerical solution corresponding 
to an X-ray flare peak energy of ~ 1 keV. In this case the vF v 
peak flux of the IC component is roughly two orders of mag- 
nitude larger than the synchrotron one; this is a good condition 
for the detection of high energy flares. In Fig.[13]we show syn- 
chrotron (solid line) and IC (dotted line) spectra ~ 200 s after 
the start of the flare, and in Fig. Q~2] we give synchrotron (1 
keV; solid line) and IC (1 GeV; dotted line) light curves, and 
compare the latter with the 1 GeV LAT sensitivity (dot-dot- 
dot-dashed line). The fluence carried by the high energy flare is 
significantly larger than that carried by the X-ray flare, and in- 
tegrating for a sufficiently long time, LAT can easily detect the 
expected high energy flare. However, in such a case the peak 
energy of the high energy flare is of the order of TeV (see Fig. 
[T3l). thus there is a temporal delay between the X-ray and the 
high energy flare in the LAT energy band. 

5. Individual cases: GRB 011121 and XRF 011030 

In this section we estimate the IC component in the framework 
of DES for XRF 01 1030 and GRB 01 1 121. In the case of XRF 



01 1030, lGalli& Piro 1 (12006) have shown that a late X-ray flare 
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Fig. 11. Synchrotron (solid line) and IC (dotted line) spectra 
for a thick shell fireball expanding in an ISM ~ 200 s after 
the start of the flare, fo ~ 500 s. The model parameters are 
£53=10.0, r =150, «=10.0, e e =0.01, e B = 10~ 3 , p = 2.5 and 
Z = 1.0. With this set of model parameters the fireball is in 
the fast cooling regime at the time of the flare occurrence, and 
the mean flare peak energy is around 10 eV. The dot-dot-dot- 
dashed line is the 5cr LAT sensitivity for an integration time of 
500 s. The vertical dot-dashed line indicates the cutoff energy 
due to pair production for the IC spectrum. 
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Fig. 13. Synchrotron (solid line) and IC (dotted line) spectra 
for a thick shell fireball expanding in an ISM ~ 200 s after 
the start of the flare, to = 500 s. The model parameters are 
£,3 = 10.0, T =150, «=10.0, e c =0.27, e B = 10~ 5 , p = 2.5 and 
Z - 1.0. With this set of model parameters the fireball is in 
the fast cooling regime at the time of the flare occurrence, and 
the mean flare peak energy is around 1 keV. The dot-dot-dot- 
dashed line is the 5<x LAT sensitivity for an integration time of 
500 s. The vertical dot-dashed line indicates the cutoff energy 
due to pair production for the IC spectrum. 
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Fig. 12. Synchrotron (solid line) and IC light curves at 100 
MeV (dotted line) for a fireball expanding in an ISM in the 
framework of DES. The model parameters are the same as in 
Fig-QTJ The origin of the time is shifted to f =500 s. With this 
set of model parameters the fireball is in the fast cooling regime 
at the time of the flare occurrence, and the mean flare peak en- 
ergy is around 10 eV. The dot-dot-dot-dashed curves shows the 
5cr 100 MeV LAT sensitivity as a function of the integration 
time. 
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Fig. 14. Synchrotron (solid line) and IC light curves at 1 GeV 
(dotted line) for a fireball expanding in an ISM in the frame- 
work of DES. The model parameters are the same as in Fig. 
IT31 The origin of the time is shifted to ?o=500 s. With this set 
of model parameters the fireball is in the fast cooling regime at 
the time of the flare occurrence, and the mean flare peak energy 
is around 1 keV. The dot-dot-dot-dashed curve shows the 5cr 1 
GeV LAT sensitivity as a function of the integration time. 



occurs in this event around 1300 s and it can be explained in the 
context of the DES scenario, both for a fireball expanding in a 
wind environment and in an ISM. However, in their computa- 



tion, Galli & Piro ( 20061) did not account for the cooling due to 
the presence of an IC component, thus now we want to assess 
whether the presence of this component affects their findings. 
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In the ISM case, iGalli & Piro I (12006ft found a broadband 
(radio-to-X-ray) solution for a jetted fireball with £53 = 0.03, 
T = 130, n = 5, s e = 0.29, e B = 8 ■ 10~ 5 , p = 2.1, z = 1 
and break time Tt — 8 • 10 5 sec. This solution implies that dur- 
ing the fast cooling regime the relative importance of IC and 
synchrotron components is Y ~ 27, thus the X-ray synchrotron 
flux is suppressed by a factor (1 + Y) ~ 30. On the contrary, 
with this set of model parameters, the optical and radio bands 
are below the cooling frequency v c at the time of observations, 
thus in these bands the synchrotron flux is not affected by the 
presence of IC emission. To increase the predicted X-ray flux 
in order to be consistent with the observations, one should in- 
crease the fireball energy £53 and/or eg, but this has the effect 
of also increasing the predicted optical and radio fluxes bring- 
ing them above the observed fluxes (see Eqs. in appendix B of 
Panaitescu & Kumar (2000)). We thus conclude that in pres- 
ence of IC emission for a fireball expanding in an ISM, we 
cannot find a broadba nd solution. 

In the wind case, IGalli & Piro I d2006l) found a broadband 
solution with £53 = 0.3, r =60, A«=0.055, e e =0.02, e B =0.001, 
z=l and p=2.1. With such parameters, in presence of an IC 
component the X-ray synchrotron flux is reduced at maximum 
by a factor ~ 3 during the fast cooling regime. At the time of the 
flare the fireball passes from the fast to the slow cooling regime, 
thus v x > v c > Vj and £ v> . c oc v~ p,2 vf~ m vl /2 oc (1 +Y)~ l . In this 
case the cooling frequency increases with time, thus around 10 6 
s v c > v x and F vx does not depend on v c . Finally, optical and 
radio emission are also not affected by the presence of an IC 
component. This means that in this case, only the X-ray light 
curve around the time of the flare occ urrence is affected by the 
presence of IC emission. As shown in IGalli & Piro (2006), in 
this phase the emission strongly depends on the initial fireball 
Lorentz factor To, thus we can find a new broadband solution 
simply by increasing this parameter from Fq = 60 to To = 70 
(see Fig. \T5t . The presence of an IC component also has the 
effect of reducing the cooling frequency v c by a factor (1 + Y) 2 , 
thus causing a delay of the time when it crosses the X-ray band. 
However, the spectral break related to the crossing of v c in the 
X-ray band is still consistent with that suggested by the two 
late time CHANDRA observations (red squares in Fig. ITSb. In 
Fig. Q3] we show the predicted synchrotron (blue line) and 1 
GeV IC (light blue line) light curves. As expected from the ES 
model, there is a good temporal correlation between the X-ray 
and GeV flare. Indeed, at the time of the X-ray flare the peak of 
IC component is around 1 GeV, and this also ensures that the 
presence of an IC component does not affect the spectrum of 
the X-ray flare. We also note that even if XRF 01 1030 is not a 
bright burst, the predicted IC emission related to its X-ray flare 
is potentially detectable by the LAT. 

Another interesting case is that of GRB 011121. The broad- 
band afterglo w data of GRB 1112 1 and its X-ray flare were 
explained bv lPiroetain d2005l) and IGalli & Piro I d2006l) by a 
thick shell fireball expanding in a wind-like medium, with ef- 
ficiencies e c =0.01 and eg =0.5. As we have showed in Sect. [2] 
IC emission can dominate over synchrotron emission only if 
e e > £b- This means that even if the X-ray flare occurring in 
GRB 011121 is very bright, a non significant contribution of 
IC emission is expected for this event. 




Fig. 15. Synchrotron (blue) and 1 GeV IC (light blue) light 
curves for XRF 01 1030 in the case of a thick shell fireball ex- 
panding in a wind density profile. The model parameters are 
£53=0.3, r =70, A«=0.055, e e =0.02, e B =0.001, and p=2.1.The 
black dot-dashed line is the LAT sensitivity at 1 GeV as a func- 
tion of the integration time. 



6. Discussion and conclusions 

In this paper we study the relative importance of synchrotron 
and IC emission in the framework of ES for a fireball expand- 
ing in an ISM and in a medium with a wind density profile. We 
first focus our attention on the "standard" case of a thin shell 
fireball, and study synchrotron and IC emission as a function 
of the model parameters that mainly determine the relative im- 
portance of these two emission mechanisms. We find that the 
importance of the IC component increases with the (fe/eg) ra- 
tio during both the fast and slow cooling regime, and that for 
a given(e e /fB) value it is greater during the fast cooling regime 
than the slow cooling regime. We then evaluate IC flux as a 
function of e e and eg and we find that for typical GRB model 
parameters during the first phases of afterglow emission, e.g. 
~500 s, if the fireball is in the fast cooling regime, IC flux in- 
tegrated over the full energy range can reach values as high as 
~ 10~ 8 erg cm~ 2 s~ l . However, at early times, the peak of the IC 
component can be at very high energies, well above the GeV 
range, depending on the position of the peak of the synchrotron 
component. We thus study the detectability of IC component 
as a function of e e and eg in the Large Area Telescope (LAT) 
above 1 GeV as a function of the integration time. We find that 
for a sufficiently long integration time, such as ~ 10000 s, if 
the fireball is still in the fast cooling regime LAT can detect IC 
afterglow emission for a significantly larger region of param- 
eters space. We develop a numerical model which permits to 
study ES synchrotron plus IC emission at times smaller than 
the fireball deceleration phase. We produce complete early-to- 
late time synchrotron and IC light curves and spectra, and com- 
pare them with the LAT sensitivity. As suggested by the stud- 
ies performed in the e^, (q/eb) parameters space, when most 
of the emission occurs during the fast cooling regime LAT can 
detect IC emission from the afterglow, and the detectability of 
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this component increases with the (e e leg) ratio. One can easily 
extrapolate these results to the case of AGILE, where the effec- 
tive area around 1 GeV is a factor of ~ 15 lower than LAT. The 
main conclusion is that we would expect detection by AGILE 
up to a redshift z ~0.25; this is consistent with the findings of 
Gou & Meszaros I J2007I) . 

The large IC fluxes that can be produced in the afterglow 
phase motivated us to apply this model to the very bright GRB 
940217, which presented high energy emission lasting from 
~ 500 s to ~ 5000 s after the burst. We find that IC emission 
from ES can account f or this delayed emission. We note also 
that IWei & Fan I (120071) recently proposed an alternative sce- 
nario, in which this delayed high energy emission is ascribed 
to IC emission arising from a decelerating blast wave in the 
presence of a late energy injection plus temporal evolution of 
the shock parameters e e and eg- 

The case of DES (thick shell fireballs) produced by a long 
lasting central engine activity is interesting, because it offers a 
possible explanation for the origin of flares with a soft spec- 
trum consistent with that of the subsequent afterglow emission, 
or that of flares showing a spectral evolution AT of the order of 
0.5-1.0. In this context X-ray flares represent the onset of af- 
terglow emission, thus the flare photons can be IC scattered by 
the afterglow electrons, producing flares in the GeV-TeV band. 
We calculate synchrotron and IC spectra and light curves, and 
find that bright GeV flares can be produced for a burst with 
fireball kinetic energy of the order of ~ 10 54 erg if the fireball 
is in the fast cooling regime at the time of the flare appear- 
ance. The analysis of a sub-sample of bursts with known red- 
shift having flares in their X-ray light curve shows that ~ 30 % 
of them satisfy this requirement (see Table|2]l, and are expected 
to present detectable high energy flares. Short GRBs typically 
have lower isotropic energies, see e.g GRB050724 in Table [2] 
thus we do not expect detectable high energy flares from this 
population of bursts. A comparison of model predictions with 
the LAT sensitivity shows that LAT is able to detect high en- 
ergy flares and IC emission arising from the afterglow up to 
~ 1 GeV if most of the emission occurs during the fast cooling 
regime, and that the detectability of the IC component increases 
with the (e c /eB) ratio. Given an X-ray flare with mean fluence 
~ 3 x 10~ 7 erg cm" 2 , if the Compton parameter Y is of the 
order of unity, the expected high energy flare can be detected 
by LAT at 100 MeV, while its fluence is below the 1 GeV LAT 
threshold. However, we have shown that for typical model pa- 
rameters Y can be larger than unity, and in this case the fluence 
of the high energy flare overcomes the 1 GeV LAT threshold. 
We also estimate, under reasonable assumptions, the values of 
flux and peak energy of the expected high energy flare from 
those of the X-ray flare. This is a very important information 
because it allows us to understand whether an observed X-ray 
flares with certain properties can produce a detectable high en- 
ergy flare. 

In the DES scenario low and high energy flares are pro- 
duced by the same electron population and source region, thus 
they are expected to be simultaneous once one considers their 
integrated broadband flux. This is an important prediction and 
represents one of the distinctive elements between the several 
models proposed to explain the origin of flares. If X-ray flares 



are produced by LIS associated with a long-lived central en- 
gine, then high energy flares can be produced through two dif- 
ferent mechanisms. One possibility is that X-ray flare photons 
are IC scattered on the electrons accelerated by the forward 
shock (External Inverse Compton, EIC). In this case, due to the 
fireball curvature and the anisotropy of the incoming photons 
in the frame of forward shock electrons, the high energy flare is 
expected to extend up to longer time scales than the X-ray flare. 



and th e flux received by an observer is lowered (see lFan et al. 



(2007) for details). The second possibility is that high energy 
flares are produced in LIS by SSC of the same electrons pro- 
ducing X-ray flares. In this case one expects a good temporal 
correlation between the X-ray and high ene rgy flare, as in th e 
framework of DES. However, as shown by iFan et al. I (120071) . 
in the framework of LIS, SSC emission is expected to peak at 
lower energy with respect to the ES model. One should also 
note that to observe the simultaneity of low and high energy 
flares, the peak energies of both flares need to fall within or 
below the observational band of the instrument. We predict 
that X-ray flares with peak energy of ~ 10 eV produce high 
energy flares with peak energy around hundred of MeV-GeV. 
Consequently, to distinguish between LIS and ES models re- 
quires both spectral and temporal information. Therefore fu- 
ture GLAST observations, and their coordination with those 
of the Swift satellite, will play an important role in narrowing 
down models attempting to explain the origin of flares, and thus 
shedding light on the physics of the central engine. 
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